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The precision of biochemical signaling is limited by randomness in the diffusive arrival of molecules 
at their targets. For proteins binding to the specific sites on the DNA and regulating transcription, 
the ability of the proteins to diffuse in one dimension by sliding along the length of the DNA, in 
addition to their diffusion in bulk solution, would seem to generate a larger target for DNA binding, 
consequently reducing the noise in the occupancy of the regulatory site. Here we show that this 
effect is largely cancelled by the enhanced temporal correlations in one dimensional diffusion. With 
realistic parameters, sliding along DNA has surprisingly little effect on the physical limits to the 
precision of transcriptional regulation. 



INTRODUCTION 



gives fewer independent samples. 



Cells constantly regulate the expression levels of their 
genes. A central motif in this regulatory process is the 
binding of transcription factor proteins to specific sites 
along the DNA. The precision of transcriptional regula- 
tion is limited, ultimately, by randomness in the arrival 
of transcription factor (TF) molecules at these sites pQ. 
But proteins can find their binding sites on DNA by two 
very different mechanisms — either by diffusing in three 
dimensions through the surrounding solution, or by bind- 
ing weakly and diffusing in one dimension along the con- 
tour of the DNA molecule. The idea of diffusion or slid- 
ing along DNA goes back (at least) to the realization that 
the lac repressor seems to bind more rapidly to its target 
site than would be allowed by three dimensional diffu- 
sion alone [2J|3]. More recently, the discussion of 3D vs. 
ID diffusion has been revitalized by theoretical analysis 
of the optimal search strategies [4], by new biochemical 
measurements [5] , and by direct physical observations of 
the sliding motion [HI 13 E] • Here we consider the impact 
of dimensionality and diffusion on the physical limits to 
the precision of transcriptional regulation. 

Any physical system which responds to the concen- 
tration of a signaling molecule will exhibit noise due to 
the random diffusion of these molecules in and out of 
its "sensitive volume" [T]; the larger the volume, the 
smaller the (fractional) noise. This leads to the intu- 
ition that ID diffusion will have a huge impact on the 
noise because it effectively increases the size of the tar- 
get to which the transcription factor is binding [9 . Our 
main result is that this intuition is wrong. The prob- 
lem is that diffusion in one dimension has a very differ- 
ent statistical structure than in three dimensions, and 
so one cannot simply say that ID sliding generates a 
larger 3D target. We show that realistic combinations 
of 3D and ID diffusion in fact have a surprisingly small 
effect on the limiting noise level — the increased size of 
the target is largely cancelled by stronger temporal cor- 
relations, which means that integrating over a fixed time 



II. CONCENTRATION FLUCTUATIONS 

Before starting on the problem of proteins binding to 
DNA, it is useful to recall some facts about concentra- 
tion fluctuations for molecules free in solution [10 . If we 
write the concentration as c(x,t), then the power spec- 
trum of fluctuations S c (\l,uj) is defined by 



(<5c(x,t)<5c(x',0) 



du> 



S' c (k,cj)e 



+ik-(x-x')-iw(t-t') 



,(1) 



where d is the dimensionality of space. For molecules 
present at the mean concentration c and diffusing freely 
with diffusion constant D, the power spectrum is given 

by 

„ , 2cDk 2 

* c(k ' w)= (z?fcY + - 2 ' (2) 

where k — |k|. This result is independent of the dimen- 
sionality d, although of course the units of concentration 
are different in different dimensions. If we integrate over 
all frequencies, which corresponds to making instanta- 
neous measurements, we find a spatial power spectrum 

g,S c (k,u;) = c. (3) 

This is spatial white noise, and embodies the fact that a 
snapshot of molecules in solution will reveal a random, 
Poisson distribution: the variance of the molecule num- 
ber in any finite volume is equal to the mean, so the 
power spectrum of concentration fluctuations is equal to 
the mean concentration, and again this is independent 
of dimensionality. 

Consider a measurement of concentration that is av- 
eraged over some small region of space, 



C(t) = J d d xT^(x)c(x,i), 



(4) 
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where W is a weighting or windowing function that de- 
fines the size of the region we are interested in. To in- 
terpret C as an average of the local variable c, we must 
have J d d xW(x) — 1. The temporal fluctuations in C 
are themselves determined by a power spectrum Sc(ui), 



where 
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The normalization of W implies that W{k = 0) = 1. If 
the size of the region is £ d , then VI 7 (k) will fall to zero for 
k much larger than A ~ l/l. Thus we can approximate 
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We are interested in the behavior of the power spec- 
trum at low frequencies, corresponding to measurements 
with long averaging times. If we try simply to set u = 0, 
then we find 



S c (w = 0) 



For d = 3, 
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This corresponds to white noise in the time domain, so 
we expect that averaging over time will reduce the noise 
variance in proportion to the averaging time. More pre- 
cisely, if we average for a time Tint, we will be sensitive 
to frequencies \u\ < l/rj n t, so we will see a variance 



((5C) 2 ) 3d 



M<l/Ti„t 

r 

M<l/T in t 

c 



dio 
2^ 
cko 
2?r 



S% d (uj = 0) 



DIt Sv 



(12) 
(13) 
(14) 



Although we have not done a full calculation including 
the effects of binding and unbinding to target sites, this 
is essentially the "noise floor" for a system which senses 
the concentration of signaling molecules by having them 
bind to a site of size £ [T] fTT] [12]. 
In contrast, for d = 1 we have 
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which is infrared divergent. To get the right low fre- 
quency behavior in one dimension we have to be a bit 
more careful. We have 
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This integral is actually finite as A 
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oo, so we can write 




Now the variance of measurements averaged over a time 
Tint becomes 
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We see that in one dimension, the variance of concen- 
tration measurements declines only as the square root 
of the measurement time. This is contrast to the three 
dimensional case, where the standard deviation of the 
noise declines as the square root of time but the vari- 
ance declines in direct proportion to the measurement 
time, as seen in Eq (14 1. 



III. BINDING WITH 3D DIFFUSION 

If transcription factors bind specific sites on the DNA 
to control the rate of gene expression, then the noise 
in the regulated gene product will contain a contribu- 
tion from the noise in the occupancy of the specific site. 
Here we want to compute this noise contribution, start- 
ing with the case where the transcription factors find 
their cognate site by 3D diffusion alone and progressing 
towards more complicated translocation strategies in the 
following sections. 

In this section we briefly review the calculations of 
Ref [TT] . Consider a cell volume with an average concen- 
tration c of TFs present, and a binding site located at 
x = 0. The transcription factor molecules can diffuse 
with a bulk diffusion constant D$, bind to the specific 
site at a rate k + c, and dissociate back into the free so- 
lution at a rate k-, according to the following set of 
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equations: 



dc 
&t 



= D 3 V 2 c-M(x-x ), 



h = fc + c(xo, t)(l — n) — k- 



(24) 
(25) 



To determine the fluctuations in the binding site occu- 
pancy n, we linearize the equations about the mean oc- 
cupancy n(t) = Sn(t) + n with n = k + c/ (k + c + fc_), and 
around the mean concentration, c(x, t) = 5c(x,t) + c, 
and write the perturbations as Fourier modes: 



Sn(t) 



duo 

W) 1 
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We are looking for the power spectrum of fluctuations in 
occupancy, which we can define by 



(28) 



If there exists a low frequency limit, S n (uj — > 0), and 
we average the noise for a time Ti nt that is longer than 
all other relevant timescalcs in the problem, we can [by 
analogy with the arguments leading to Eq (14)] use it 
to calculate the observable variance by computing er 2 = 
S n (u> -> 0) /rint- 
Having introduced the necessary notation, we can 
compute the fluctuations in concentration at the binding 
site, x = 0, from Eq pi) : 



(5c(xq, uj) = iujSn 



d 3 k 
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This integral is ultraviolet divergent, and must be cut 
off at A = — , where a is the binding site size; keeping 
leading terms in iu> one obtains: 



<5c(x , uj — > 0) = 



2TrD 3 a 



(30) 



The linearization of Eq ( 25 1 yields 



Sti - 
- iujSh = h fc+(l - n)<5c(x , uj) + k-fi/36F, (31) 

where r^ 1 = k + c + fc_ is the time scale of occupancy 
fluctuations, /3 — 1/ksT, and SF is the fluctuation in 
the free energy difference between the bound and empty 
states of the site; it is the 'force' that is thcrmodynam- 
ically conjugate to the 'displacement' Sn. After substi- 
tuting the local concentration fluctuations, Eq (30 1, into 



Eq (31), one can finally use the fluctuation-dissipation 



theorem to compute the power spectrum of the noise in 
occupancy: 
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The result of Ref [TT] was that, in the case of 3D diffu- 
sion to the binding site, the low- frequency limit of the 
occupancy noise power spectrum is: 
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The first term corresponds to the binomial switching 
fluctuations as the occupancy of the specific site changes 
between full and empty; this term depends on the micro- 
scopic details of the TF-DNA interaction (here embodied 
by the off-rate constant The second term is caused 
by the fluctuating diffusive flux to the binding site, and 
provides a lower bound on the noise, independent of de- 
tails. Expressed in terms of the equivalent concentration 



fluctuations, S c (u>) = S n (u>) 
for the second term is 
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the noise spectrum 



S c (uj —>■ 0) = c/irD^a, 



with the associated fractional variance 
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where, to facilitate comparison with results that follow, 
we explicitly pull out the leading numerical factor and 
1 /(effective binding site size). This is the noise lower 



bound, consistent with the result of Eq ( 14 1 for I = a, 
that depends solely on the rate of 3D diffusion D3, the 
binding site size a, and the amount of time averaging 
that the system performs, T- lnt . 



IV. BINDING WITH ID DIFFUSION 

Consider the case where the diffusion in bulk does not 
occur, and the TF is only free to slide along the length of 
the DNA. Let £,(x, t) denote the ID concentration along 
the DNA contour, which we imagine to be stretched in 
the x direction, and let D% be the corresponding diffusion 
constant. Equations ( 24|25 1 remain unchanged but for 
the new notation, i.e. c(x, t) — > £(a;,t),Z?3 — * D\, V 2 — > 
<9 2 ; we will use q to denote the Fourier variable conju- 
gate to x. We again need to calculate the concentration 
fluctuations at the binding site (cf . Eq ( 29 1 ) : 



S^xq, uj) = iujSh 



dq 



(27T) 



(36) 



This integral diverges as q — > if we set w = 0, but in 
contrast to the discussion above there is no problem as 
q — > 00. The full result is that: 
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Note that unlike Eq (30 1, this result does not contain 



the binding site dimensions. Using Eq (37 1 we compute 
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the diffusive contribution to the noise power spectrum 
(cf . Eq ((34| , but no limit ui 



0): 



(38) 



The noise variance is obtained by integrating the noise 
power spectrum, as in the arguments leading to Eq ( 23 1 , 
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V. COMBINED ID AND 3D DIFFUSION 

In the discussion so far we have been missing a param- 
eter that would interpolate between two qualitatively dif- 



ferent noise regimes: the 3D result of Eq (35 1 and the 
ID result of Eq ( 40 1 . If the rate for dissociation from the 



sliding state into the bulk is increased, the ID case must 
approach the 3D result; the residence time will grow ever 
shorter and will, at some scale, break the strong corre- 
lations reflective of the ID sliding mode with its high 
probability for returning to the origin. 



We see that the noise variance declines as 1 / y/T mt , con- 
sistent with the simpler calculation of concentration fluc- 
tuations in Eq ( 23 1 . This again is in contrast to our intu- 
ition that variances should decline as 1 / r;„t , which is the 
result for binding coupled to 3D diffusion. The differ- 
ence between ID and 3D can be understood by realizing 
that in ID diffusion, a particle leaving the binding site 
at the origin has a large probability (in fact, probabil- 
ity 1) of returning back to the origin. A receptor trying 
to estimate the local concentration of TFs will therefore 
be unable to collect samples that are truly independent, 
and the variance in the measurements will consequently 
decrease at a rate that is slower than expected from av- 
eraging over independent measurements. 

Suppose that the DNA with the sliding TFs is embed- 
ded into the cytoplasm, where the bulk concentration of 
TFs is c; the TFs can jump onto the DNA at a rate (per 
unit length) k+c, where we expect k+ ~ D3. The TF will 
stay on the strand for a residence time, r r = k_, where 
k_ is the rate for dissociating from the DNA. There will 
be an equilibrium between the ID and 3D concentra- 
tions, c = £/T r «+ « ^/r r D3. With these identifications, 
the noise variance can be rewritten in a form similar to 



A. Perfect mixing in bulk solution 

We first solve a simplified case where the TF diffuses 
on the DNA, but if it dissociates into the bulk, it mixes 
very quickly and there is no correlation between the 
point of dissociation and subsequent re-association to the 
DNA strand. We expect that the divergence observed in 
the pure ID case, Eq (36 1, will now get regularized. The 
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the result in Eq ( 35 ) 
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Here, t he eff ective binding site size is a c g = 
\/DiT r \A" r /Ti n t. Naively, one could think that the ef- 
fective site size would be equal to the average length 
that the TF explores during ID diffusion on the DNA, 
b = y 'D\T T , and therefore a c g ~ 6, but that would be 
wrong: the effective site size depends on r int to compen- 
sate for the highly correlated fluctuations in ID diffu- 
sion. Because Tj nt ^> r r [13], the effective site size will 
be much smaller than b. The apparent reduction in the 
noise expected in the naive picture - because the size 
of the binding site, a, is replaced by presumably much 
larger length explored by ID diffusion, b - must in real- 
ity be traded off against the longer required integration 
time. As a result, it is not immediately clear whether 
ID diffusion decreases the noise lower bound compared 
to 3D case. 



binding site 
occupancy 
n{t) at xo 



FIG. 1: The transcription factors can either be free in so- 
lution at concentration c, or they can enter a region on the 
DNA where they diffuse by sliding. The ID concentration is 
denoted by £(x,t). The specific binding site on the DNA is 
at location xo = 0; k+£(xo) and fc_ are the on- and off-rates 
for transitions between the ID sliding state and the bound 
state on the specific site; k+c is the rate per unit length for 
transition onto the DNA from free solution, and k_ = t^ 1 is 
the rate for dissociating from ID solution into the 3D solu- 
tion. The effective radius of the DNA molecule is R and the 
"sliding length," or the average distance along the contour 
covered in a single ID random walk before dissociation, is 
denoted by b. 
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dynamical equations for this case are: 
dn 



dt 

dt 

dC 
~dt 



= k+£(xo, t)(l — n) — k^n 

d 2 i 



(41) 



Dl %^ -h8{x-x ) - K-£+ k + C (42) 



k_ J dx £,{x,t) — k + LC 



(43) 



Here, C(t) is the bulk concentration that is perfectly 
mixed and therefore does not have any spatial depen- 
dence; TFs in the bulk are enclosed into a box with 
a side of length L and the DNA inside is an infinitely 
thin line of length L; and «_ and n + are rates for 
transitioning between ID concentration £ and the bulk 
C. When linearized and written out in terms of their 
Fourier components, the equations couple to the zero 
mode, 5l;(q = 0,u>), differently than to non-zero modes: 

(-w + l/T c )6n = k+(l - n)5£(x ) + k-n[38F(4A) 
(-iw + Dig 2 )5| 9 = iuSh- K-S£ q + K+L5C (45) 
{-iuj + k + L)SC = K-<5£ g =o (46) 



We use Eq (45 1 at q — together with Eq (46 1 to express 



SC and substitute it back into Eq (45 1 for non-zero q; it 



is then straightforward to write S^(xo) as: 
dq iljY,(uj) Sn 



(2n) -iu> + Diq 2 + k_ 



where 
E(cj) = 1 



K-K+L 



(— iuj + «;_)(— iuj + n+L) — K-K+L 



(47) 
(48) 
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Note that as u> — > 0, S(w) will diverge; before computing 
the power spectrum of fluctuations S n using Eq ( 32 1 we 



J 



therefore need to expand Eq (47 1 to the leading power 
of iw. On the other hand, the momentum integral in Eq 



We describe the system by the following set of equations 

dn 

— = k + €{x ,t){l - n) - k-i 
dt 

dt(x,t) _. d 2 £(x,t) dn 



(47) now converges due to the k_ term in the denomi- 



nator. 

After some algebra, the effective power spectrum in 
the concentration fluctuations is: 



S((u -» 0) 



2 v /L»ik_ 



(l + 7) ; 



(50) 



where 7 = k + L/k- is the ratio of on- and off-rates for 
going into ID solution. The 7-dependent expression in 
parenthesis has a minimum of 7/8 at 70 = 3, and there- 
fore a noise floor exists that cannot be eliminated by 
some convenient combination of the rate parameters. 

Rewritten in terms of the "equivalent concentration" 
c = C/L 3 = k_^/k + L 3 w k_^/D 3 , the concentration 



noise floor is (cf. Eq (40 1) 
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- = — x x . (51) 

CJ 16 D 3 CT int _ y/DiT r 

This is the result we might naively have expected, 
namely that sliding along the DNA creates a larger tar- 
get, of size a e g ~ b. But we get this result only because 
we have assumed that the TF can repeatedly dissoci- 
ate into a well-mixed bulk solution, instantly losing the 
memory of the point along the DNA from which it dis- 
sociated. 



B. Diffusion in bulk solution 

We will now solve the coupled 3D- ID diffusion prob- 
lem without assuming that in the bulk solution the tran- 
scription factors mix perfectly. The relevant quantities 
are schematized in Fig[T] 



dt 

dc(x, R, t) 
dt 



£>i- 



d 2 a 



dt 



S(x — xq) + k + c(x, R = 0, t) — K-!;(x, t), 



D 3 V 2 c(x, R, t) - [k + c(x, R = 0, t) - n_£(x, t)] S(y)S(z). 
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We have again assumed that DNA is stretched along the x-axis and that it is an infinitely thin molecule. £ is a 
function of only one variable, x, while c is a function of coordinate x and radial coordinates R that are perpendicular 
to the x direction. 

We linearize and Fourier transform the equations as follows: 
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where again q is the spatial Fourier variable conjugate 
to x and kj^ is conjugate to R. Note that 5l; q is function 
only of q, while Sc q .k ± is a 3D Fourier transform of the 
bulk concentration fluctuations that depends on both q 
and kj^. Finally, Sc(q,R — 0,tu) are Fourier modes of 
concentration fluctuations along the x direction, evalu- 
ated at the location of the DNA strand, R = 0. 



We first compute Sc q from Eq (57 1 
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Here, k 2 = q 2 — iuj/D 3 , 



and A = -| is the ultraviolet 
cutoff at the radial size of the DNA molecule. We can 
substitute Sc q from Eq ( 60 ) into Eq ( 56 1 to obtain the 
expression for 5£(x ): 



5£(xo) — ituSn 



dq 
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Compared to the pure ID result, Eq (36 1, Eq (61 1 now 
contains a new term k__F^ 1 (<7) in the denominator; for 
k_ = this term vanishes and the result reverts to the 
ID case as expected. F(q) depends on the momentum 
q, whereas in the perfect mixing case, Eq (47 1, it was 
simply equal to 1. The integrand of Eq (61) is still di- 
vergent for oj — as q — > 0, but the integral nevertheless 
converges. Assuming a non-zero the limit u> — > 
therefore exists, and we can integrate: 



5i(x o ,u^O)= lU ^I(a,0), 
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The noise power spectrum for occupancy follows in close 
analogy to Eq ((Ml: 
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The effective spectrum of local concentration fluctua- 
tions corresponding to the second term is (cf. Eqs (34 



38l[50l) 
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The associated noise variance is: 
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With the above definition of a, the equilibrium be- 
tween ID and 3D concentrations turns out to be c = 
Air D 3 (xt y . The noise variance in concentration can be 
rewritten as: 
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(70) 



The result looks similar to the pure 3D diffusion case 
in Eq ( 35 1 : the noise due to concentration fluctuations 



has its length scale a (receptor size in pure 3D case) re- 
placed with R (effective DNA cross-section in the com- 
bined ID / 3D case) , and the noise contribution term gets 
multiplied by the pref actor that we examine next. 

We expect the parameter a = n + /4irD 3 to be close to 
1 for diffusion limited approach to DNA. Imagine that 
the area that TF attempts to hit in order to bind non- 
specifically to the DNA is a cylindrical segment of DNA 
with radius R and length b. If we were treating the 




FIG. 2: The relative decrease in noise magnitude, compared 
to the pure 3D diffusion model, as a function of parameters 
a and /3. Three values for a are shown, spanning two orders 
of magnitude (a — 1 solid, a — 0.1 dotted, a = 0.01 dashed); 
f3 covers the relevant range if typical ID diffusion length is 
as expected from search time optimality arguments (order 
hundred base pairs). 
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cylindrical DNA segment as a sphere of radius 6, then 
Smoluchowski on-rate limit ~ AnD^b would apply |20j . 
In the first approximation, the on rate per unit length, 
k + , would then be AttD 3 and a = 1. Exact deriva- 
tions for Smoluchowski limit in simple geometries were 
reviewed in Ref [2T]; for elongated objects (prolate el- 
lipsoid with long semiaxis b and short semiaxis R), the 
diffusion limited rate is given by AnD^b/ ln(2b/R), and 
therefore a ~ ln _1 (26/i?). a is hence less then one, and 
has a weak logarithmic dependence on the ratio of the 
scale of the persistence length and the effective radius of 
the DNA. 

The parameter (3 is approximately the square of the 
ratio between the cross-section of the ID cylinder (the 
"target" that 3D diffusion has to hit) and the average 
"sliding length" b along the DNA. R must be of order of 
several nanometers; while b is, at DNA stacking length 
of a = 0.3 nm per base-pair and 100 bp average diffusion 
length jU |T3], around b = 10 — 100 nm. It is therefore 
not unreasonable to assume that the factor /3 could be 
as low as (3 ~ 10 -3 — 10~ 2 . The corresponding decrease 
in noise relative to pure 3D diffusion is shown in Fig [2] 
Note that as (3 — > oo one should recover the pure 3D 
diffusion result. Looking back at expression for I{a,j3) 



in Eq (64 1, as (3 increases, the t 2 term in the denominator 
of the integrand is becoming irrelevant. If we neglect it 
completely, the integral is solvable analytically: 



I(a,/3-> oo) = 0- 



A 



(71) 



where A^ = - is the cut-off at the binding site size along 
the DNA, and A x /A = R/a. Inserting the integral in the 
large (3 limit into the noise result, Eq (70 1, we see that 
(3 cancels and we get: 



1 / R \ 1 1 

7T H X X — 

2n V a a J D 3 cn nt R 



> 



1 



1 



1 



2ira D 3 cT int a 



(72) 
(73) 



which is essentially the 3D result, Eq (35 1 



VI. DISCUSSION 

While the role of ID diffusion of transcription factors 
along the DNA has been recognized as a possible expla- 
nation for the observed fast DNA-TF association rates, 
the impact of this additional mode of TF translocation 
on the noise in binding site occupancy remained unex- 
plored. The question is important for two reasons: first, 
the diffusive contribution to the noise in gene expression 
must fundamentally limit the precision of transcriptional 
regulation; and second, there is an appealing argument 
that ID diffusion could drastically increase the target 
area on the DNA that TFs have to find and correspond- 
ingly lower the limiting diffusive noise. Here we show 



that this intuitive argument is wrong — while there might 
be some reduction in the noise if the bulk diffusion is sup- 
plemented by ID sliding, this reduction is not expected 
to be significant. 

Much has been said about possible TF translocation 
strategies on the DNA such as three-dimensional vol- 
ume exchange, local dissociation-reassociation reactions 
(hopping), sliding along the DNA and intersegmental 
transfers mediated by DNA- looping and about 
the ways in which these mechanisms influence the tar- 
get search times [U [331 1211 While our model does 
not examine all of the proposed mechanisms and makes 
approximations about the geometry of the TF and the 
DNA, we are ultimately only interested in what happens 
locally around the specific site as opposed to computing 
global properties such as target search time statistics. As 
a result, while rare but long range excursions might be 
important in search arguments, they should not signifi- 
cantly affect the noise as long as they are not the dom- 
inant form of all bulk transfers; if they are, one would 
have the case discussed in Section V.A of perfect mix- 
ing upon dissociation. Moreover, the conclusions pre- 
sented here reflect the basic differences in the diffusive 
processes in one and three dimensions, in particular the 
high probability of returning to the origin in ID diffu- 
sion, and are therefore not a consequence of the detailed 
assumptions about the TF-DNA interaction, as is clear 
from the "back of the envelope" arguments presented in 
the Section II. 



The diffusive contribution to the total noise in gene 
expression has only recently been recognized as signifi- 
cant. Metzler discusses the concept of "interaction vol- 
ume" around the regulatory site on the DNA and studies 
the probability that transcription factors will enter this 
volume in the case of the A-phage infecting Escherichia 
coli bacterium, concluding that the spatial fluctuations 
can be important in genetic circuits |16j . Holcman et 
al study the master-diffusion equation in the context 
of signaling molecules binding to and unbinding from 
ion channels [17] • The analysis of Bicoid/Hunchback 
system in the development of the fruit fly Drosophila 
melanogaster shows that the measured signatures in 
the noise in Hunchback expression are consistent with 
diffusive fluctuations in Hunchback's regulator, Bicoid 
[20 HH]. Van Zon and coworkers [TS] study the diffu- 
sional component of the noise in Green's function reac- 
tion dynamics (GFRD) stochastic simulations and con- 
clude by noting that ID sliding along the DNA could 
have important effects on the noise power spectra. All 
these results point to the basic physics of diffusion as 
setting a limit to precision of a fundamental biological 
process, transcriptional regulation. Our results indicate 
that this fundamental limit is not easily evaded by the 
(still largely unknown) complexities of protein motion 
along the DNA molecule. 
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